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Abstract. We study the collapse towards the gravitational radius of a macroscopic 
spherical thick shell surrounding an inner massive core. This overall electrically neutral 
macroshell is composed by many nested 5-like massive microshells which can bear non- 
zero electric charge, and a possibly non-zero cosmological constant is also included. The 
dynamics of the shells is described by means of Israel's (Lanczos) junction conditions 
for singular hypersurfaces and, adopting a Hartree (mean field) approach, an effective 
Hamiltonian for the motion of each microshell is derived which allows to check the 
stability of the matter composing the macroshell. We end by briefly commenting on 
the quantum effects which may arise from the extension of our classical treatment to 
the semiclassical level. 



1. Introduction 

The dynamics of collapsing bodies is a subject which attracted much attention amongst 
physicists since the very first formulation of General Relativity. However, the general 
problem of gravitational collapse is extremely complicated and simplified models have 
therefore been constructed which allow one to describe the dynamics at least for the 
simplest configurations (for a review of analytical solutions, see Ref. |T]). In this 
perspective, a lot of effort has been dedicated to the study of the collapse of thin 
spherical shells (see, e.g. Refs. 013111), for they represent a quite general and simple 
model and can be regarded as the starting point for the analysis of more complicated 
systems. In recent years, attention has been focused on quantum effects that can arise 
from strong gravitational fields characterizing the final stages of the collapse and several 
approaches have been attempted (for a very limited list, see, e.g. Refs. (SHEIIZI)- 

In particular, a formalism was introduced in Ref. jH] which allows one to study the 
semiclassical behaviour of the shell by adopting a "minisuperspace" approach jH] to its 
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dynamics very much in the same spirit as that previously employed for the Oppenheimer- 
Snyder model in Ref. [TUj. Yet, many aspects of the collapse still need to be investigated 
at the classical level, which may alter the quantum (semiclassical) picture in a substantial 
way. For example, the issue of whether the radial degree of freedom of the shell can 
be described by a quantum wave function strictly depends on the fact that the matter 
composing the shell be confined within a very small thickness (i.e. of the order of the 
Compton wavelength of the elementary particles composing the shell) around the shell 
mean radius, so that the assumption of a coherent state is reliable. We refer to this issue 
as the "stability" problem for thick shells. In fact, a realistic description of the collapse 
should involve finite quantities (that is, a shell with finite thickness) whereas, already 
for the classical dynamics, the very use of standard junction equations |3J requires 
the validity of the thin shell limit (i.e. the shell's thickness be much smaller than its 
gravitational radius). 

In this paper we address the classical dynamics of a collapsing spherical shell (which 
we call macroshell) of proper mass M with an inner massive core of mass M c and study 
the conditions under which the matter composing the shell remains confined around its 
mean radius. In order to verify this property, we discretize the continuous distribution of 
matter and consider the macroshell as if composed by a large number N of homogeneous 
nested spherical £-like microshells jllj . Exploiting the spherical symmetry, we can 
adopt a Schwarzschild-like coordinate system and the motion of the microshells can be 
described by the usual "areal" radial coordinates r, (i = 1,...,N), ordered so that 
fi < r 2 < ...rjy. Therefore, the thickness of the macroshell is given by 5 = (rjy — J"i) and 
its mean radius is R = (rjv + 5"i)/2. Each microshell will be assumed to have the same 
proper mass = m and a (possibly zero) electric charge Q{. We will limit our study 
outside the gravitational radius of the system, Rjj = 2G M s /c 2 , where M s is the total 
ADM mass, so that the radius of the innermost microshell T\ will always be greater than 
Rh- 

To summarise, the model will be constructed assuming the following conditions: 

(i) M c 3> M. The mass of the core is much bigger than the macroshell's, so that the 
total mass-energy of the system is approximately M s ~ M c ; 

(ii) m <C Nm = M. Each microshell is much lighter than the macroshell; 

(iii) (R >) ri 3> 5. The inner radius of the macroshell (and therefore the radius of 
each microshell) is much larger than its thickness. This implies that the macroshell 
thickness also remains much smaller than Rjj and this condition will limit the 
temporal range of our analysis; 

TV 

(iv) >J Qi = 0. The macroshell is overall electrically neutral; 



dR 

R(0) > R H and — 



0. The collapse starts in the weak field regime. One can 

t=o 

therefore take t as the Schwarzschild time, although it will then be more convenient 
to use the (micro)shell proper times (for a detailed comparison, see the Appendix B 
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in Ref. |TT|). 

As mentioned above, the classical equations of motion for a thin shell in General 
Relativity can be obtained from the standard Israel's junction equations. In the present 
paper, we shall adopt a minisuperspace approach in which the embedding metrics 
(interior and exterior with respect to the shell) are held fixed and chosen to be particular 
solutions of the (bulk) Einstein equations, and derive the shell's dynamical equations 
from a variational principle applied to an effective action. This (equivalent) treatment 
will allow us to take into account the matter internal degrees of freedom and can be 
extended to the semiclassical level along the lines of Refs. jHlHIl- In fact, the equations 
thus obtained are equivalent to the junction conditions, with the exception of one more 
equation which has no counterpart as a junction condition and represents a secondary 
constraint for the time variation of the gravitational mass of the shell. The relevant 
dynamical equation is a first order ordinary differential equation for the shell radius 
which contains ADM parameters (those defining the inner and outer metrics). In 
order to determine the shell trajectory, one can fix the initial shell radius, but the 
initial velocity is then uniquely determined if the ADM parameters were also fixed. 
Alternatively (and more sensibly), one can set the initial radius and velocity and adjust 
the ADM parameters correspondingly. 

The task of extending the above analysis to a system of (micro) shells would be 
extremely involved in either cases. If we choose to set the initial radii and velocities, 
then the ADM parameters will have to be reset each time two shells cross. Alternatively, 
if we set the ADM parameters once and for all, it is possible that the system runs into 
some spurious (cusp or "superficial") singularity when two shells cross [T| [13] . In light of 
a future extension to the semiclassical level, we shall therefore find it more convenient 
to derive an effective Hamiltonian for the motion of a single microshell in the mean 
gravitational field generated by the core and all the other microshells treated as a 
whole. This will allow us to obtain a one-particle effective potential which describes 
the tidal forces the microshell is subject to, and, in particular, we will see that the 
binding potential varies with time and confines the microshells within an ever decreasing 
thickness (at least sufficiently away from the gravitational radius of the system). 

We will work in units in which the speed of light c = 1 but will explicitly display 
the gravitational (Newton) constant G. 

2. Single shell dynamics 

We begin by reviewing the case of one 5-like shell discontinuity. The action for a generic 
matter-gravity system is given by 

S = Seh + Sm + Sc , (2.1) 

where Seh is the bulk Einstein-Hilbert action, Sm the matter action and Sc accounts 
for boundary contributions. 
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Figure 1. A schematic view of the relevant space-time. The surface at r = Rb serves 
to bind the collapsing system inside a "box" of finite spatial extension and the limit 
Rb — ► oo can be eventually taken. 



In order to describe the motion of a shell, let us consider a space-time volume VL 
with boundary dQ = U^E^, in which every represents a smooth three-dimensional 
submanifold (see Fig. HJ. The action takes the more explicit form (see Refs. [HUH] and, 
for a thorough discussion, Ref. [TK] ) 

S = 1 [ d 4 x y/^g 71 + / d 4 x C m 

+ 5^77 / d3x + ^ I d " x smh-\ V ) , (2.2) 

O 7T G JU fc E fc 7TU Jn fe Sj. 

where g is the determinant of the four-metric g^ u , TZ the four-dimensional curvature 
scalar, K, the extrinsic curvature of the 3-dimensional hypersurface whose metric 
has determinant ( 3 'g and r] is the scalar product of the unit normals to two (non- 
smoothly) intersecting boundaries and whose volume element is denoted by 
|( 2 )^|V 2 . Boundary terms have been included so as to obtain the equations of motion 
by varying the metric inside Q with vanishing variations on the border dfl E3 ■ 

In our case, we must also consider the three-dimensional surface S swept by 
the time-like sphere of area AttR 2 which splits the space-time Q into two regions: 
Qi n , internal to the shell, and f2 out , external to the shell (where radiation can be 
included [OH)- Of course, there are an initial hypersurface Si and a final hypersurface 
Sf with possible discontinuities at the shell trajectory's end-points S ; /f R E. Adopting 
a minisuperspace approach, the metrics inside the two regions fim/out are chosen to be 
spherically symmetric solutions of the Einstein field equations and can be written in the 
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form 



ds u 



-Adt 2 + A~ 1 dr 2 + r 2 dn 2 



(2.3) 



in which r is the usual areal coordinate and t is the time variable t\ < t < tf. The most 
general solution of this kind (without external radiation) is the Reissner-Nordstrom- 
(anti-)de Sitter metric with 



.4, 



2 CM 



in/out 



in / out 



+ 



Qin/out 



± 



A 



in/out ^2 



(2.4) 



r 3 
where Mi n / out is the ADM mass (proper mass plus gravitational energy) enclosed in the 
sphere of radius r, Qi n / ou t represents the electric charge inside the sphere of radius r, 
A-in/out the cosmological constant (the signs before the last term account for both the 
de Sitter and the anti-de Sitter metrics). 

In Ref. [El, it was shown that, upon choosing a suitable ADM foliation [12J of the 
space-time, the shell's trajectory can be expressed as r = R(r), where r is an arbitrary 
time variable on the shell's world volume E, and the three metric on E is given by 

-N 2 dr 2 + R 2 dQ 2 



dsl 



[2.5) 



where N = N(r) is the lapse function of the shell. On then integrating over space 
coordinates, the reduced effective action (|2.2j) for the shell can be written as 



S eS (N,R,M Q 



1 

G 



dr 



dr 



pR-RRt&nh 



-i 



AnNR 2 p- 



M 



2A 



out 



R 



2A 



out 



drL e{i (N,R,M out ) 



[2.6) 



where (3 = y AN 2 + R 2 , the subscripts "in" and "out" now denote quantities calculated 
respectively at r = R — e and r = R + e in the limit e — > + , p = p(R, M out ) is the 
mass-energy density of the shell, and we have introduced the notation [F]™ t = F in — F out . 

We point out that the action (|2.fij) allows for the time variation of the ADM mass 
M out which is a canonical variable of the system, and this corresponds to a radiating 
shell. In this case, the region Q ont would be further divided into two regions, one 
enclosing all the radiation and one (O^ in Fig. [TJ devoid of both matter and radiation. 
In the following Sections, the ADM mass calculated on the external surface of the shell 
is constant and therefore the last two terms in the action will be dynamically irrelevant. 
However, we will give the more general expressions for the radiating case below. 

Upon varying the action with respect to the canonical variables N, R and M out one 
obtains the Euler-Lagrange equations of motion for the shell in the form 

££ G ff _ dL cS _ d_ ( dLes" 
dr 



6X 



dX 







dX 



(2.7) 
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where X = (N,R,M ont ) and SX(t[) = SX(t{) = 0. Introducing the canonical momenta 
conjugated to the variables N, R, M out one obtains 



Pr = 



dN 
dL ea _ 
dR 
dL eS 



+ 



M n 



out 



2A 



out 



P N = ^ = (2.8) 

(2.9) 

P M = -^5E_ = — ^— . (2.10) 
dM out 2A out 1 ; 

We see from the first equation that N is a Lagrange multiplier, corresponding to the 

effective action being invariant under time reparameterization on the shell. Varying the 

effective action with respect to N will hence yield a primary constraint instead of a real 

equation of motion. Eqs. (|2.7|) thus read 



5L 



off 



5N 
5R 



R 

GiV 
1 

G 



IP] 



m 

out 



(3 + 



SL, 



eft 



-AttN 
AttNR 



P 

5R 

2 Sp 



- 4nR 2 p = 

•■ N 2 M 
RR 



(2.11) 



RNR\ 



R 



Po 



N 



R 



+ 



out 



R 



- out 







SM out (l-2GM out /R) 
At this point, one can introduce the Hamiltonian of the system 

P R R + P M M out -L eS = 47rNR z p-^ \f3\ 



/3 out l-2M aat /R 

(2.12) 

(2.13) 



H 



(2.14) 



A-kNR 2 P GL ^ Jout , 

and it becomes clear that Eq. (|2.11j) represents the (secondary) Hamiltonian constraint 
in the form 

H 

(2.15) 



N 







whereas Eq. ()2.12|) can be written as the time preservation of the constraint ()2.15|) . 



d 

d7 



o . 



(2.16) 



Making use of the invariance under time reparameterization on the shell, it is convenient 
to choose the proper time gauge iV = 1, so that the dynamical evolution of the shell is 
completely determined by the equation 

- 4 7T R 2 p 



H 







(2.17) 



and t is the proper time on the shell. If we consider a shell with constant ADM mass 
M out , it seems obvious to identify the physical shell energy with its proper mass and set 
4tcR 2 p = M. With this substitution, Eq. (j2~TTj) reads 

R ( r~. • „ i~. ■ „\ 

(2-18) 



M 



G 



R 2 - \ A 



R 2 



where we have used the explicit expression for (3 and R 
derivative with respect to the shell's proper time r. 



dR/dr denotes the total 
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Figure 2. Evolution of the macroshell's velocity with the radius for M s = M, M c = 
and A = 1 (in units with G = 1). 



2.1. Dynamical equation and initial conditions 

On squaring twice the above Eq. 1)2.18)1 and rearranging we have 



R2 = T^PTT^ - A>ut) 2 - \ (An + A out ) + ^J- . (2.19) 



It is then clear that the shell's trajectory will be uniquely specified by the initial 
condition R{r — 0) = Rq if the metric functions A in / out have fixed parameters. One is 
therefore given two options: 

a) fix the metric parameters and the initial radius -R(O) = Rq; 

b) set the initial radius -R(O) = Ro and velocity -R(O) = Ro, and adjust the metric 

parameters accordingly. 

For example, let us consider a macroshell collapsing in a Schwarzschild-de Sitter 
background, with metrics 

2GM C A „ 

A in = l + - r 2 

r 6 

(2.20) 

2GM S A 9 

A ont = 1 ~ + -r 2 , 

r 6 

where M c is the ADM mass of the core and M s the total ADM mass of the system ^f. 
Substituting Eqs. (l2~27H) into Eq. (j2~lHl) . we obtain 

• 2 ( M s -M c \ 2 G(M S + M C ) G 2 M 2 A 
R = \~^I~) ~ l + R + ^^"3^ ' (2 ' 21) 

which gives the time evolution of the shell's velocity in terms of its radius R(R(r)). 
Clearly, (as Fig. El also shows) that equation describes the collapse of a shell starting 

Note that M s — M c > is not equal to the shell's proper mass M in general. In fact, it can be smaller 
than M for bound orbits since it also contains the shell's negative gravitational "potential energy" . 
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from its maximum radius Rq = i? max (M, M s , M c , A) with zero initial velocity Rq = 0. 
In order to be able to set the initial velocity independently of the initial radius, one can 
relax one of the metric paramaters, e.g., the total ADM mass M s , which will therefore 
be determined as M s = M S (R , R ; A, M, M c ) by Eq. (l2~2TJ) at r = 0. 



2.2. Generalised dynamics 



As previously noted, one cannot set the initial radius and velocity if the metric 
paramaters have already been fixed. This, for example, rules out the possibility of 
setting M s = M c + M, A = and have the shell start collapsing with zero velocity at 
finite radius. 

Another possibility is to consider the presence of extra (non-gravitational) forces 
which certainly come into play in the description of astrophysical objects, and allow 
for equilibrium configurations of the collapsing matter at finite radius (for example, the 
radiation pressure for main-sequence stars and the degeneracy pressure of the electrons 
and, at later stages, of the neutrons for compact stars). A way to model that is to 
introduce an ad hoc force which acts upon the system before the collapse begins and 
manages to keep the radius R fixed (for r < 0). Such a force will be switched off at 
t = and the shell then collapses under the pull of just the gravitational force. One 
can study this (generalized) dynamics of the shell by adding a constant Vq to the shell 
effective Hamiltonian, 



H, 



cff 



M 

T 

M 

T 



R 2 + V{R) + V 



R 1 



where 



fo 



M 



M s - M c 
M 



[V(R ) - V 



G (M s + M c ) G 2 M 2 



R 



AR 2 



,(2.22) 



M, — Mr 



1 - 



G (M s + M c ) G 2 M 2 A 



+ 



-rR, 



M J ' R ' 4i2g 3"°' (2 ' 23) 

is now an arbitrary constant. We can now set the initial conditions R{0) = Rq and 
-R(O) = for a given Rq. Let us finally remark that the Hamiltonian constraint is still 
given by 

H eS = , (2.24) 

and the analysis then proceeds straightforwardly. In fact, Eq. ()2.22|) for a non-radiating 
shell is still compatible with the constraint (|2.16|) . since fo and all the metric parameters 
are constant in such a case and the dynamical system remains well defined. 



3. The macroshell model 



In our model, the macroshell is collapsing in a black-hole background and the total 
ADM mass of the system is constant. In order to study the stability of the matter 



Classical dynamics and stability of collapsing thick shells of matter 



9 



M 



U t. 



M, 



Figure 3. Internal structure of the macroshell. We have denoted with M < = M- ln and 
M> = M out the ADM mass calculated respectively on the internal and on the external 
surface of the microshell of radial coordinate x. 



composing the shell, its internal structure must be taken into account and therefore 
the macroshell is viewed as if composed by many nested microshells. We expect the 
matter distribution inside the macroshell to alter locally the black-hole background and 
we study the motion of a single microshell in the mean gravitational field generated by 
the remaining (N — 1) microshells. 

In order to do this, let us single out a microshell of radius x and proper mass 
m = M/N (see FigEJ). When x ^ Ti, we can make use of Eq. f)2.17j) with 4nx 2 p = m, 
substitute the expression for (3, square twice and rearrange so as to obtain 

where x = dx/dr denotes the total derivative with respect to the shell's proper time 
r. Then, on multiplying Eq. (J2.19|) by m/2, we obtain the equation of motion for a 
microshell in the form of an effective Hamiltonian constraint, induced by the canonical 
Hamiltonian constraint ()2.17|) . namely 

H eS = -x 2 + V eS (x) = , (3.2) 

and comparison with Eq. (|2.19|) immediately yields 

VW*) = = {A< + A>) - £^ (A* - A>)« - ^ . (3.3) 

In order to proceed, it would now be necessary to substitute the explicit expressions 
for A < and A>, which requires the knowledge of the microshells distribution inside the 
macroshell. As we discussed in Section |2~TI there are two equivalent ways of determining 
the evolution of x = x{t), one which allows one to set both x(0) = x and x(0) = x 
(with metric parameters correspondingly determined), and one in which the metric 
parameters are fixed and just x(0) = Xq can be set. However, the former option is 
presently further complicated by the presence of more than one shell. In fact, whenever 
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Figure 4. Equivalent configuration in the thin-shell limit: (1) x < X; (2) x > X. 



two shells cross (i.e. at the time r = r c when x(t c ) = for some % — 1, . . . , N — 1), the 
metric parameters must be re-evaluated using x(t c ) and x(r c ) as new initial conditions + . 
Since r c in turn depends on the trajectory before the crossing, this also means that the 
system keeps memory of its entire evolution. For this reason, and because we want to 
be able to extend the treatment to the semiclassical level in the future, we find it more 
convenient to pursue the other option and fix the metric parameters once and for all. 
This may imply that the generalised dynamics of Section 12.21 will be used in order to 
allow us to impose arbitrary initial velocity (e.g. in Section HJ. We shall however not 
display the arbitrary constant fo in order to keep the expressions simpler. 

We then choose the metric parameters according to the assumptions as 
outlined in the Introduction. First of all, since m <C M <C M B , we shall take 
M> — M< ~ m, which amounts to x(0) ~ and is consistent with the choice of 
having the microshells confined within a thickness 5. With this assumption both M< 
and M> become functions of the microshells distribution {(r«, Mj), i = 1, . . . , N}, but, 
since it is unknown a priori, the evaluation of the potential remains involved. A reliable 
approximation is therefore to adopt a mean field approach and study the motion of the 
microshell when x < r\ or x > r\ + 5; in the thin shell limit for the macroshell, all the 
remaining (N — 1) microshells can be viewed as forming a shell of proper mass (M — m) 
and, if we denote its average radial coordinate by X, we can obtain an equation of 
motion analogous to Eq. ([2.190 . In agreement with our gauge choice, we will use the 
proper times of the (two) shells to foliate the space-time and therefore consider the 
dynamical variables x and X as functions of the same variable r. We remark that, if 
the microshell bears an electric charge Q, then the shell of mass (M — m) must have a 
charge — Q so that the total macroshell of mass M has zero charge. 

One has therefore two distinct situations for x < X and x > X (see cases (1) and 

+ We do not intend to discuss the issue of shell crossings in any detail here. What happens to two 
shells when they collide would in fact require the implementation of specific matter properties (such as 
the detailed treatment of non-gravitational interactions) . 
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(2) in Fig. BJ and two distinct effective potentials: 

For x < X, the region Rh < r < x has a Schwarzschild-de Sitter metric with ADM 
mass (M s — M) , the region x < r < X has a Reissner-Nordstrom-de Sitter metric with 
ADM mass (M s — M + m) and electric charge Q, while the outer region r > X has 
again Schwarzschild-de Sitter metric with ADM mass M s . We therefore get 



m 
~2 



G 



x 



(2 M s - 2 M + m) + 



Q 2 



G 2 m 2 



Gmx Ax 2 



+ 



2x< 



Q 2 



2G 2 m 2 



A 2 
+ 3* 



(3.4) 



M 



Hi 



- — (2M s -M + m) 
X 



Q 2 



G 2 (M 



m) 



G(M — m)X 



4X 2 



+ 



_Q1 

2A 2 



U^x 2 

2G 2 (M-m) 2 3 



(3.5) 



where Eq. f!3.4|) represents the effective potential for the microshell of proper mass m 
and radius x, while Eq. ()3.5j) is the effective potential for the macroshell of proper mass 
(M — m) and radius X. 

The same procedure can be followed for x > X. In this case the region r < X has 
a Schwarzschild-de Sitter metric with ADM mass (M s — M), the region X < r < x has 
a Reissner-Nordstrom-de Sitter metric with ADM mass (M s — m) and electric charge 
— Q, while the outer region r > x has Schwarzschild-de Sitter metric with ADM mass 
M s . This yields 

Q 2 G 2 m 2 



cff 



X 



m 
~2 



-- (2 M s 

x 

Q 2 



m 



Gmx Ax 2 
Q 2 \ A 



2 a; 2 



2G 2 m 2 



3 X 



(3.6) 



KffPO 



M-m 



G 
X 



(2 M s — M — m) + 



Q 2 



G 2 (M - m) 



G(M-m)X 



4 A 2 



+ 



2X 2 



1 




(3.7) 



2G 2 (M- 

Since Eq. (13 .2j) must hold for both the microshell and the macroshell, for each one of 
the two cases presented above we can write a total effective Hamiltonian for the system 
which satisfies the constraint 



H%(x,X)=H* {[ (x) + HX(X) 



m 



M 



rn 



in 



ir + 



M 



m 



X 2 + V: s (x) + V*(X) 
X 2 + V^(x,X) = . 



(3.f 



2 2 

The dynamics has now been reduced to that of a two-body system, so that we can 
separate the motion of the center of mass from the relative motion by introducing the 
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center of mass and relative coordinates given by 
f = x — X 

(3.9) 

„ m u ,„ 

R =M X + X ' 

M m 

where \l = m (M — m)/M is the reduced mass of the system. The Hamiltonian 
constraint ()3.8|) then takes the form 

H ir) {f,R) = |f 2 + yi? 2 + ^ (T) (r,i?) = , (3.10) 

where the superscript r is to recall the time dependence of the Hamiltonian through the 
variables f and R. 

Since one component of the system (the macroshell) is much more massive than the 
other one (the microshell) we have, to lowest order in m/M, \x ~ m and we may thus 
consider the motion of the microshell in the center of mass system to be described by 
the relative coordinate f in much the same way as one treats the motion of the electron 
in the hydrogen atom. Indeed, since the macroshell is freely falling in the black-hole 
background, its (mean) radius R and proper time r define a locally inertial frame. 
We would thus like to subtract the center of mass motion and obtain a one-particle 
Hamiltonian for the microshell. In so doing, we obtain an effective Hamiltonian for the 
whole macroshell of proper mass M and radius R by making use of Eqs. (j3.2j) and ()3.3|) . 
In this case we just have two regions, both with Schwarzschild-de Sitter metric, r < R 
with ADM mass (M s - M) and r > R with ADM mass M s , 

H^(R) = fR 2 + vlP(R) 



• 2 2G / M\ M 2 G 2 A , 

R 2 I M s — + - R 2 

R V 2 J AR 2 3 



(3.11) 



M 

~ Y 

Finally, we obtain the one-particle Hamiltonian constraint for the microshell as 

f#>(f) = F«(f, R) - HthR) ^ f f 2 + V£\?) = , (3.12) 

where the potential may contain a term f so as to ensure that f(0) ~ 8/2 and 
f(0) ~ 0, according to the discussion of the generalised dynamics in Section l2~2l Further, 
the variable R in the potential is now to be considered as a (time-dependent) 
parameter. We note here that this approximation is justified if the typical times of 
evolution of the macroshell radius R are much larger than the period of a typical 
oscillation of the microshell around R. 

The explicit expression for the potential V^\f) is extremely complicated and we 
deem pointless to write it explicitly. Instead, in the next Section, we will plot its 
behaviour in situations of particular interest. 



4. Shell stability 



In order to study the stability of the shell we suppose that for r = all the matter is 
localized within a thickness 5 Rh R and the shell has negligible (initial) velocity 
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Figure 5. Effective potential Vm (r)/m for a macroshell thickness 5 = 6 • 10 9 Rh 
and macroshell radius R = 10 Rh- The metric parameters are M c — 3M©, M ~ 
8.5 • 10" 9 M c ~ 2.5 • 1CT 8 M Q , m ~ 4 • 1CT 9 M ~ 1 • 1(T 16 M©, M s = M c + M, and 
A = Q = 0, (M Q ~ 2 • 10 33 g is the solar mass and R H ^ 9 • 10 4 cm). 

-R(O) ~ 0. This corresponds to a microshell confined in the region |f| < 8 with maximum 
kinetic energy 

y ^ Ki T) (r = ±5) - V#>(f = 0) , (4.1) 

provided the value V^ T '(±5) is less than the (possible) maximum of Vj^\ 

In Fig. we plot the effective potential for a typical case (see the caption for the 
values of the parameters) at the (relatively) large value of R = 10 Rh and macroshell 
thickness 5^6- 10~ 9 Rh (we have chosen arbitrary constants so that Vffl(±5) = 0). 
From the graph, it appears that the microshell is classically bound around the "centre 
of mass" of the macroshell, in agreement with previous results to first order in the 
small \f\/R expansion jllj . This means that if the microshell is initially placed close 
enough to the macroshell, the mutual common gravitational attraction will keep them 
together. However, for |f| 3> 8, the potential bends down and starts decreasing. This 
effect cannot be seen in a perturbative expansion for small \f\/R but is encoded in the 
exact potential. One can understand such an effect on considering that if the microshell 
started with a radius sufficiently larger or smaller than the macroshell's, the mutual 
attraction would not be sufficient to overcome the difference in collapsing velocities and 
their "distance" |f| would increase along the collapse. We can note in passing that, 
as a consequence of the form of the effective potential, one expects that the quantum- 
mechanical probability for the microshell to tunnel out the well (i.e., from |f| < 8 to 
|f| 3> 8) is not zero and (approximately) symmetrical for the cases of a microshell falling 
inwards or outwards. Strictly speaking, bound states of microshells would therefore be 
just quantum-mechanically meta-stable. 
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Figure 6. Effective potential Vm /m for the same parameters as in Fig. and 
R=10R H (solid line), R = 8 R H (dashed line), R = 6R H (dotted line), and R ! = 4 R H 
(dash-dotted line). 

During the collapse, the velocity of the macroshell R increases while its radius R 
shrinks and the effective potential for the microshell thus changes. This behaviour 
is easily understood as the increasing of tidal forces between different parts of the 
macroshell. One might therefore wonder how the time evolution of the binding potential 
affects the macroshell's thickness and if, eventually, the microshells can escape thus 
spreading out the system. From the plots in Fig. |H1 we can see that the potential 
yij) remains symmetric around the centre but appears to become steeper and steeper 
as the shell approaches the gravitational radius. In fact, to a given maximum kinetic 
energy (|4.1|) . there correspond smaller and smaller values of 5 for decreasing values of 
R (moreover, it is easy to see that S /R actually decreases, since 5 decreases faster than 
R). On the other hand, one may also note that the "width" of the forbidden region 
decreases (whereas the maxima of the potential increase) and it is thus possible that, at 
the quantum level, the meta-stable states have different lifetimes. A definite conclusion 
on this point would need a more complete semiclassical analysis beyond the scope of 
the present work. 

For completeness, we plot the evolving potential when opposite electric charges are 
carried by the microshell and macroshell in Fig. [7[ and when a cosmological constant is 
included in Fig. |H1 

5. Conclusions 

We have examined the collapse of a macroscopic self-gravitating spherical shell of 
matter in the presence of a massive core in the framework of General Relativity. The 
(macro) shell is assumed to be formed by a large number of thin (micro) shells initially 
bundled within a very short thickness. The system is then shown to remain classically 
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Figure 7. Effective potential Vm /m for the same parameters as in Fig. [3] except 
that Q = 10 29 e (e being the electron's charge) and R = 10 Rh (solid line), R = 8 Rh 
(dashed line), R = 6 Rh (dotted line), and R — 4Rh (dash-dotted line). 




Figure 8. Effective potential Vm /m for the same parameters as in Fig. except 
that A = 10 -56 cm -2 and R = 10 R H (solid line), R = 8 R H (dashed line), R = 6 R H 
(dotted line), and R — ARh (dash-dotted line). 

stable, since the effective potential which drives the microshells becomes more and 
more binding as the macroshell approaches the gravitational radius of the system (but 
remains sufficiently away from it so that our approximations R — Rh ^> 5 and \R\ < 1 
hold). Further, including a (reasonably valued) cosmological constant does not alter 
this behaviour, nor does the presence of electric charges carried by the microshells. 

A semiclassical analysis may however change the above conclusion, since the bound 
states of microshells are just meta-stable and there may be significant changes in their 
lifetimes (or other, non-adiabatic, effects) as the macroshell collapses. In Refs. |11|IT%]. it 
was shown that the microshells can indeed be described by means of localized quantum 
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mechanical wave functions, however in the linear approximation (first order in \ f\/R) for 
the effective potential. We have now found that, for sufficiently large values of |f| (but 
still much smaller than R), the complete effective potential reaches maxima and then 
decreases, thus suggesting a non- vanishing probability for the microshells to tunnel out. 
It was also found in Refs. [HI HE] that, upon coupling the microshell wave functions 
to a radiation field, the shell spontaneously emit during the collapse as a non-adiabatic 
effect. A further analysis of such semiclassical effect in light of the present results is 
therefore in order, which could also affect the usual thermodynamical description [T7j . 

Let us conclude by mentioning that inclusion of (or a possible relation with) the 
Hawking radiation [19J as the shell approaches the gravitational radius still remains to 
be investigated. 
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